UNCLASSIFIED 


AD  NUMBER 

AD911536 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  only;  Test  and  Evaluation;  MAY 
1973.  Other  requests  shall  be  referred  to 
Department  of  the  Air  Force,  Directorate 
of  Operational  Requirements  and 
Development  Plans,  Attn:  RDQSD, 
Washington,  DC  20330. 

AUTHORITY 

USAF  ltr,  14  Jan  1974 


THIS  PAGE  IS  UNCLASSIFIED 


/&>  f// si  6 


j 


This  Division  Note  reports  research  sponsored  by  the  Directo¬ 
rate  of  Operational  Requirements  and  Development  Plans, 
DCS; Research  and  Development,  Hq  USAF,  under  Contract 
F44620-69-C-00M.  It  does  not  necessarily  represent  the 
opinion  of  anyone  other  than  the  author  and  the  Division 
Manager.  It  may  he  withdrawn  at  any  time. 


anskr  is  an  independent,  nonprofit  research  corporation  established  in  1958. 
The  Corporation  contributes  to  the  security  and  public  welfare  of  the  United 
States  by  performing  systems  analysis,  operations  research,  and  development 
planning  studies.  Its  principal  mission  has  been  and  is  to  satisfy  the  con¬ 
tinuing  need  of  the  United  States  Air  Force  for  objective  and  timely  research 
studies  of  advanced  systems,  anser  also  performs  some  research  for  other 
Government  agencies  and  nonprofit  organizations.  In  addition,  anser  con¬ 
ducts  a  small  program  of  self-sponsored  research  in  the  public  interest. 


DEFENSE  DIVISION  NOTE 


DDN  73-2 


QUICK  EVALUATION 
FOR  AIR 


COMPUTER  MODELS 
DEFENSE 


May  1973 

Victor  V.  Gogolak 


Approved  by 
R.  S.  Timm,  Manager 


- - —  DISTRIBUTION  LIMITED  — - - 

to  U.S.  Oownmnii  ooty;  t**t  *i\d  *v*lu*tion; 


Miy  1973.  Oih*  r*qo»tt*  lor  thii  mutt  lx 

«r*rr«(J  to  Mtj  USAF  JftPytDQSD). 


ABSTRACT 


This  Note  discusses  three  time-sharing  computer  programs 
that  are  quick  tools  for  air  defense  engagement  analysis. 

The  programs  are  based  on  a  Monte  Carlo  method.  The  model 
MONTYX  evaluates  air  defense  interceptor  configurations  by 
determining  their  effectiveness  against  a  bomber  cell  in 
the  terminal  phase  of  an  intercept.  The  model  DLMNTY  expands 
the  capability  of  MONTYX  to  a  situation  where  two  bomber 
types  are  contained  within  the  bomber  cell.  The  results  are 
bomber  probabilities  of  survival  for  various  interceptor -to- 
bomber  ratios.  The  model  TOOTH  involves  penetration  by  four 
types  of  bombers — two  carrying  gravity  weapons  and  two  carry¬ 
ing  standoff  weapons.  The  defense  is  provided  by  two  types 
of  interceptors  in  two  waves.  In  the  first  wave,  only  raid 
count  by  cell  is  provided.  The  second  wave  assumes  that  * 
bomber  count  has  been  made  and  standoff  weapons  have  been 
launched.  The  results  indicate  the  numbers  of  bombers  of 
each  type  surviving  each  wave  of  defense  as  well  as  the  num¬ 
ber  of  standoff  weapons  launched  and  those  which  survive  the 
second  wav©  of  defense. 
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I .  INTRODUCTION 


Models  developed  for  air  defense  engagement  analysis 
range  from  the  exponential  interceptor /bomber  (I/B)  formulas 
to  massive  and  expensive  simulations  of  the  strategic  air  war. 
Some  of  these  examine  enemy  penetrations  of  U.S.  air  defenses 
and  others  examine  U.S.  bomber  penetration. 

The  models  presented  here  have  been  in  use  for  over 
five  years  and  have  not  grown  in  complexity  or  sophistica¬ 
tion.  That  is  not  their  purpose.  They  are  quick  tools 
which  fill  the  gap  between  the  formulas  and  the  war  game 
simulations.  As  analytic  tools,  they  all  require  a  certain 
amount  of  homework  to  establish  the  inputs.  In  this  way, 
some  fairly  detailed  analyses  can  be  made,  with  the  programs 
taking  the  burden  of  the  calculations. 

The  programs  are  written  in  XBASIC  and  can  be  run  on 
any  time-sharing  computer  system  having  an  XBASIC  compiler. 

All  programs  are  based  on  a  Monte  Carlo  method.  This  allows 
programs  to  handle  such  problems  as  imperfect  commitment 
with  only  a  few  decision  statements.  Closed  form  solutions 
are  not  used. 

Applications  have  varied  from  parametric  analysis  of 
the  probability  of  detection  and  conversion  (P  )  and  probabil- 
ity  of  kill  (P  )  for  interceptors  to  assessment  of  equivalent 
megatonage  (EMT)  delivered  on  targets  in  massive  exchanges. 

A  fair  amount  of  "art"  is  involved  in  setting  up  the  inputs 
to  allow  the  analyst  close  control  of  the  games. 
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II.  SINGLE  WAVE  BOMBER  INTERCEPTOR  MODEL— MONTYX 


A.  Purpose 

The  purpose  of  the  model  called  MONTYX  is  to  evaluate 
air  defense  interceptor  configurations  (aircraft,  missiles, 
and  fire  control  systems)  by  determining  their  effectiveness 
against  a  bomber  cell  in  the  terminal  phase  of  an  intercept. 

B.  Assumptions 

The  bombers  and  interceptors  involved  in  the  interaction 
are  assumed  to  have  the  same  configuration;  i.e.,  the  bombers 
represent  identical  targets  and  the  interceptors  have  the 
same  values  for  probabilities  of  detection  and  conversion 
and  of  kill.  The  simulation  begins  after  initial  detection 
has  taken  place;  the  interceptors  are  assigned  to  bombers 
on  a  random  basis.  The  group  of  interceptors  may  make  any 
specified  number  of  passes  (P)  against  the  bomber  cell  and 
are  assumed  to  carry  at  least  P  weapons.  An  interceptor 
may  fire  one  weapon  per  pass  but  retains  its  weapon  if  it 
has  been  assigned  to  a  bomber  already  killed.  The  model  is 
run  for  various  intercept or /bomber  ratios.  The  game  may  be 
played  many  times  to  improve  statistics.  The  number  of  times 
the  game  is  played  is  specified  in  the  input  and  the  results 
are  averaged. 

C.  Program  Description 
1.  Input 

When  the  program  is  run,  the  following  numbers  are 
requested  from  the  user: 
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•  Number  of  passes  (P) 

•  Probabilities  of  detection  and  conversion  and  of 
Hill  for  each  pass:  "1"  P(DC),  P(K);  "2"  P(DC),  P(K) . 

•  interceptor-to-bomber  ratios — low,  high,  and  step 
size  II,  12,  19,  where  II  is  the  lowest  I/B  ratio 
and  12  is  the  highest  ratio  and  19  is  the  increment 
which  steps  the  game  from  II  to  12.  For  example, 

(II,  12,  19)  =  (1,  2,  .5)  would  indicate  that  the 
game  would  be  run  for  I/B  =1,  1.5,  2. 

•  Loop  size  (Ml) — which  indicates  the  number  of  times 
the  game  is  to  be  played;  the  final  results  are  the 
averages  of  the  Ml  games. 

2,  Simulation 


The  simulation  begins  after  the  inputs  are  specified  by 
the  user.  The  model  assigns  a  number  to  each  bomber.  On 
each  pass,  each  interceptor  is  randomly  assigned  to  a  bomber; 
that  is,  an  interceptor  randomly  chooses  a  bomber  to  attack. 
This  bomber  has  not  been  killed  on  a  previous  pass  but  could 
have  been  killed  earlier  in  the  same  pass  by  another  intercep¬ 
tor.  The  bomber  is  considered  killed  or  not  hilled  by  a 
random  choice  based  on  the  probabilities  of  detection  and 
conversion  and  of  kill  involved.  At  the  end,  the  total  num¬ 
ber  of  bombers  surviving  is  tallied.  After  the  specified 
number  of  passes,  some  interceptors  (those  which  did  not 
convert  on  a  target  during  one  or  more  of  the  passes)  will 
still  have  weapons.  These  are  left  unused.  The  model  pro¬ 
vides  a  number  of  runs  based  on  the  i nterceptor-to-bomber 
ratios  (I/B)  and  prints  out  the  bomber  probability  of  survival 
(Pg)  corresponding  to  the  I/B  ratio. 
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3.  Output 


if 

£ 

Ei 


The  output  consists  of  the  initial  conditions  of  the 
game  (specified  by  the  user)  and  the  bomber  (P„)  for  the 
various  I/B  ratios. 

4.  Coding 

•  Program  List  for  MONTYX 

•  Instructions  for  MONTYX 

•  Variable  List  for  MONTYX 
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PROGRAM  LIST  FOR  MONTYX 


100  DIM  FC50)#B<»0> 

110  PRINT  -NUMBER  0F  PASSES?**! 

120  INPUT  P 
130  FOR  PI  si  T0  P 

140  PRINT  "PASS**! PI  I  **P  (DC  >  #P  <K  >  ?'*  1 
150  INPUT  DCP1 )#L(P1 ) 

160  NEXT  PI 

1  TO  PRINT  -1/3  RANGE  (LOW# WISH# STEP)?” I 

180  INPUT  1 1 *12*19 

182  PRINT  "LOOP  SIZE?”! 

184  INPUT  Ml 

198  FOR  1*10*11  TO  10*12  STEP  10*19 
200  S*0 

230  FOR  J*1  TO  Ml 
240  FOR  K»I  TO  10 
250  B(KJ«! 

260  NEXT  K 

270  81*10 

280  FOR  K*  t  TO  P 

290  FOR  K1 *1  TO  I 

300  F (KI )*0 

310  NEXT  KI 

320  FOR  K 1  *  1  TO  I 

330  IF  RND(O)»0<K)  THEN  370 

340  K2*INT(t0*RND(0)«t ) 

3S0  IF  9CK2)*0  THEN  340 

360  F€K|)*K2 

310  NEXT  Kl 

380  FOR  K 1  * »  TO  I 

390  IF  F«K|)*0  THEN  460 

400  LET  X*F(K1> 

410  IF  8CX)«Q  THEN  460 

420  IF  RND(0»LOO  THEN  460 

430  9<X>*0 

440  8t«ai~t 

450  IF  B 1*0  THEN  510 

460  NEXT  KI 

470  NEXT  K 

490  S*S*8I 

510  NEXT  J 

520  PRINT  i/tO*S/<IO*Ftl> 

530  NEXT  I 
540  PRINT 

550  PRINT  -ANOTHER  GAME?”! 

560  INPUT  GS 

570  IF  Gt*-YES-  THEN  110 
1000  END 
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INSTRUCTIONS  FOR  MONTYX 

100  Dimensions  the  bomber  and  interceptor  arrays. 

110-184  Input  section;  values  of  variables  are  requested  and 
entered  from  the  terminal. 

110-120  Number  of  passes. 

130-160  Loop  indexed  by  number  of  passes  for  entering 
probabilities  of  detection  and  conversion  and 
of  hill  for  each  pass. 

170-180  I/B  range. 

182-184  Loop  size  for  iteration  of  games  for  each  I/B. 
198-530  Outer  loop  indexed  for  I/B. 

2C0  Sum  of  surviving  bombers  is  set  to  0. 

230-510  Iteration  loop  for  better  statistics. 

240-260  Bomber  array  is  set  to  1  (all  bombers  alive) . 

270  Total  bombers  in  a  game  initialized  to  10. 

280-470  PASS  LOOP. 

290-310  Interceptors  initialized  to  0. 

320-370  Interceptor  assignment  loop, 

330  Check  if  detection  and  conversion  occur, 

340  Pick  a  random  number,  1-10,  and  assign 
an  interceptor  to  that  bomber . 

350  If  bombc.  is  dead,  pick  another  number. 

360  Set  interceptor  array  to  bomber  number. 

370  NEXT  ASSIGNMENT. 

380-460  ATTACK  LOOP.  Indox  by  interceptor. 

390  If  interceptor  unassigned,  next  interceptor 

410  If  bomber  already  killed,  next  interceptor. 

420  Check  if  kill  occurs.  (If  not,  next 
interceptor.) 
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Set  bomber  to  0. 

Decrease  total  bombers  by  1. 


440 

450  If  all  bombers  in  that  game  are  dead, 
begin  next  iteration. 

460  NEXT  ATTACK 

470  NEXT  PASS 

490  Add  surviving  bombers  to  total  survivors. 

510  NEXT  ITERATION 

520  Prints  out  l/B  and  probability  of  survival. 

530  NEXT  I/B 

550-570  Check  for  another  game. 


VARIABLE  LIST  FOR  MONTYX 


F(*) 
B  {*) 
PI 
K 
K1 

I 
J 

II 
12 
19 
Ml 
P 

D  ( *) 
L  ( *) 
S 


Interceptor  Array 
Bomber  Array 
Loop  Index 
Loop  Index 
Loop  Index 
Loop  Index 

Loop  Index 
Lowest  l/B 
Highest  I/B 
Step  size  for  I/B 
Loop  size 
Number  of  passes 

Probability  of  Detection  and  Conversion  Array 

Probability  of  Kill  Array 

Sum  of  surviving  bombers  within  an  l/B  run 


5.  Sample  Runs  of  MONTYX 

NUMBER  0F  PASSES??  1 

PASS  1  P(DC)#P(K>??  .5#. 5 

I/B  RANGE  CL0W#HIGH#STEP  >??  .2*2#. 2 

LOOP  SIZE??  100 


•  2 

•  98 

•  4 

•  913 

•  6 

.878 

•  8 

.837 

i  * 

.788 

1.2 

.737 

1  .4 

.728 

1.6 

.684 

1  .8 

.629 

2. 

.611 

ANOTHER  GAME??  YES 

NUMBER  CF  PASSES??  1 

PASS  I  PCDO#PCKm  .7**5 

I/B  RANGE  (L0W#HIGH#STEP>  ??  .2#2#.2 

LOOP  SIZE??  100 


•  2 

•  968 

•  4 

.879 

.  6 

.815 

•  8 

.786 

1. 

.729 

1.2 

•  664 

1  .4 

.614 

1.6 

.61 

1  *8 

.543 

2. 

.506 

ANOTHER  GAME??  YES 

NUMBER  0 F  PASSES??  1 

PASS  I  PCDC>#PCK>??  .S*.7 

I/B  RANGE  CLOW# HIGH# STEP ) ??  .2#2#.2 

L00P  SIZE??  100 


.2 

.965 

.4 

*908 

.6 

•  835 

*8 

.782 

1. 

.724 

1.2 

.681 

1.4 

.623 

i  .6 

.575 

1.8 

•  568 

2. 

.532 

ANOTHER  GAME??  YES 
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NUMBER  0F 
PASS  1 
PASS  £ 
PASS  3 
l/B  RANGE 
L00P  SIZE 
.2 
.4 
.6 
*8 
I. 


PASSES??  3 
P <QC ) #P (K  >  ??  .7**7 
P<OC)#P<K>??  .6#. 6 
P<OC>*P(KJ??  •  Sfb 
<L0W»HIGH*STEP) ??  .2*8#.2 
7?  100 

« 894 
.679 
.497 
.348 
.202 


1.2 


.127 


1  .4 
1.6 
1  *8 
2. 


5.50000E-02 

2.20000E-02 

1.0000QE-02 

4.00000E-Q3 
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SAMPLE  RUSS  USING  KC83TYP — Continued 
(Merge  MOSTYZ,  MONT YX ,  MONTYP) 
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III.  SINGLE-WAVE  DUAL  BOMBER -INTERCEPTOR  MODEL— DLMNTY 


A.  Purpose 

The  purpose  of  DLMNTY  is  to  expand  the  capability  of 
MONTYX  to  consider  a  situation  wherein  two  bomber  types  are 
contained  within  the  bomber  cell. 

B .  Program  Description 

1.  Simulation-Input -Output 

The  assumptions  and  operation  of  DLMNTY  are  identical 
to  those  of  MONTYX;  the  only  differences  are  in  the  input 
and  output. 

The  following  numbers  are  requested  as  input  from  the 

users 

•  P  «  number  of  passes. 

•  W  <=  bomber  scaling  factor.  For  example,  W  »  2 
indicates  there  are  twice  as  many  bombers  of  type 
A  as  of  type  B. 

•  Loop  size  M,  which  indicates  the  number  of  times  the 
game  is  to  be  played. 

•  Probabilities  of  detection  and  conversion  and  of  kill 
for  each  pass  against  each  bomber  type. 

•  Interceptor-to-bomber  ratio — low,  high,  and  step 
size.  17,  19,  and  18,  where  17  is  the  lowest  I/B 
ratio,  19  is  the  highest  ratio,  and  18  is  the  incre¬ 
ment  which  steps  the  game  from  17  to  19. 

The  output  consists  of  the  initial  conditions  of  the  games 

(specified  by  the  user)  and  the  probabilities  of  survival 

(P  )  for  the  total  bomber  cell  as  well  as  the  P  for  each 
b  S 

bomber  type. 


PROGRAM  LIST  FOR  DLMNTY 


- 

y 

£  • 

k 

>;  “ 
f 


{ 


100  DIM  F (2Q0)*B(50> 

110  DIM  AC 50) 

120  PRINT  ’’WHAT  ARE  PASSES*  BOMBER  SCALE  FACTOR*  AND  LOOP  SIZE?**) 
130  INPUT  P  *  W  #  M 1 
140  IF  P* 0  THEN  1010 

150  PRINT  "NEW  PROBABILITIES  ARRAY?(YES-l *N0-O)"l 
160  INPUT  Y 
170  IF  Y=0  THEN  220 
180  FOR  Pl=l  TO  P 

190  PRINT  "PASS" I  PI J"P<DC)*P<K>  C30MBER  A) , PCDC ) *P(K )  (BOMBER  9)"l 
200  INPUT  CCP1 )*M(P1 ) *DCP1 >*LCP1 ) 

210  NEXT  PI 

220  PRINT  "1/8  RANGE  <L0W*HIGH*AND  STEP)?"! 

230  INPUT  17,19,18 

240  PRINT  TA3C 1 0  )  J "I /3"l TAB (25) J "PROBABILITY  OF  SURVIVAL" 

250  PRINT  TAB (29 ) I "ALL  BOMBERS  TYPE  A  TYPE  3" 


310 

FOR  1=17*10 

*(W*1)  TO  I9*lO*(W+l)  STEP  I8*10*CW*1> 

320 

S*0 

330 

T«0 

340 

FOR  J=1 

TO 

Ml 

350 

FOR  K*1 

TO 

10 

360 

3 (K )«l 

370 

NEXT  K 

380 

FOR  K»1 

TO 

!Q*W 

390 

A  (K  )  a  1 

400 

NEXT  K 

410 

31  =  10 

420 

AlaW*10 

430 

FOR  K=l 

TO 

P 

440 

FOR  K l *  I 

TO 

I 

450 

F (K 1 )*0 

460 

NEXT  HI 

470 

FOR  K 1  a  1 

TO 

I 

480  IF  RN0(0)>S/(!*W)THEN  545 
485  IF  81*0  THEN  600 
490  IF  PND (0 )>D (K )  THEN  600 
510  K2® INT C 1 0*RND (0 )*l ) 

520  IF  BCK2  )*Q  THEN  510 
530  FCK1 )“K2 
540  GO  TO  600 
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PROGRAM  list— Continued 


545  IF  A1 =0  THEN  600 
550  IF  RN0(0>>C<K)  THEN  600 
570  K3=INTCW*10*RNDC0>*1  > 

580  IF  A  <K3 )*0  THEN  570 

590  FCKl )=1 00+K3 

600  NEXT  K 1 

610  FOR  Kiel  TO  I 

630  IF  F<KI>«0  THEN  770 

630  IF  F<K1  )>100  THEN  700 

640  LET  X=F(K1) 

650  IF  B<X)»0  THEN  770 
660  IF  RND  (0)*»L<K )  THEN  770 
670  B  < X >  s0 
680  B! «8 I *1 
690  GO  T0  750 
700  LET  Z*F<KI >-100 
710  IF  A (Z 5 *Q  THEN  770 
720  IF  RN0<0) >M(K)  THEN  77 0 
730  LET  A<Z)*0 
7  40  A 1  °A1  <*  1 
750  IF  B1<>0  THEN  770 
760  IF  Aleo  THEN  850 
770  NEXT  Kt 
780  NEXT  K 
800  S«S*91 
830  T*T*A1 
850  NEXT  J 
860  01*10*O*W> 

870  Q2«M1*QI 
880  93*I0<*M1 
890  04«W*03 

900  PRINT#IOOO#I/01#<S^T>/02*T/Q4#S/Q3 

950  NEXT  I 
960  PRINT 
970  PRINT 
980  PRINT 
Q90  00  TO  120 

1000  FMT  F9«2#X24#F7»3#X6#F7*3#X6#F7*3 

1010  ENO 
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INSTRUCTIONS  FOR  DLMNTY 


100-110  Dimension  F,  B  and  A. 

120-230  Input  Section. 

120-130  Input  number  of  passes  for  interceptors,  bomber 
scale  factor  (10  times  V  is  number  of  bombers 
of  type  A  in  each  cell)  and  loop  size, 

140  If  0  was  entered  for  P,  then  game  is  over  and 

program  goes  to  END. 

150-170  Program  asks  if  a  new  probabilities  array  is  to 
be  set  up. 

180-210  If  YES,  then  a  loop  is  set  up,  indexed  by  the 
number  of  passes,  to  input  the  probabilities 
of  detection  and  conversion  and  of  kill  for 
each  of  the  two  types  of  bombers. 

240-250  Output  headers, 

310-950  Outer  loop,  indexed  by  I/B. 

320-330  Sum  of  surviving  bombers  of  each  type  set  to  0. 

340-850  Iteration  loop  indexed  up  to  the  specified  loop 
size. 

350-400  Both  bomber  arrays  B  and  A  are  set  to  1  (all 
bombers  alive) . 

410-420  B  bomber  total/cell  set  to  10?  A  bomber  total 
cell  set  to  10*w. 

430-780  PASS  loop. 

440-460  Interceptors  initialized  to  zero. 

470-600  Interceptor  assignment  loop. 

480  Pick  a  random  number  and  assign  the 

interceptor  to  the  B  or  A  type  of 
bomoer  based  on  the  ratio  of  the 
bombers. 
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INSTRUCTIONS  FOR  DLMNTY— Continued 


485  If  type  B  all  killed,  next  interceptor. 

490  Check  if  a  type  B  bomber  is  detected. 

510-530  Randomly  assign  the  interceptor  to  a 
live  type  B  bomber. 


545  If  type  A  all  killed,  next  interceptor. 

550  Check  if  a  type  A  bomber  is  detected. 

570-590  Randomly  assign  the  interceptor  to  a 

live  type  A  bomber  {ICO  is  added  so  that 
the  assignment  can  later  check  for 
type  A  or  B. 


610-770 

620 

630 

650 

660 

670 

680 

700 

710 

720 

730 

740 

750-760 

770 


Attack  loop. 

If  interceptor  was  not  assigned,  return 
interceptor. 

If  assigned  to  type  A,  go  to  700  (i.e.,  if 
assignment  number  is  over  100)  . 

If  bomber  already  killed,  next  interceptor. 

Check  for  kill  (Monte  Carlo) . 

Set  bomber  =0. 

Decrease  type  B  bomber  total  by  1. 

Subtract  100  to  discover  the  type  A  bomber 
assigned  to  the  interceptor. 

If  all  killed,  next  interceptor. 

Check  for  kill. 

Set  bomber  -  0. 

Decrease  type  A  total  by  1. 

Check  if  all  bombers  of  each  typo  killed. 

If  YES,  exit  to  next  iteration. 

End  ci  attack  loop. 


780  End  of  pass  loop. 

800-830  Add  surviving  bombers  to  respective  totals,  by 
type. 
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INSTRUCTIONS  FOR  DV.MNTY — Continued 


950 

990 

1010 


850  Next  iteration. 

860-890  Calculation  of  output  scaling  factors  determined 
by  W  and  Ml. 

900  Print  overall  l/B,  total  probability  of  survival 

P  for  each  type, 
s 

Next  I/B 

Go  to  beginning  to  check  for  a  new  game. 

END. 


VARIABLE  LIST  FOR  DLMNTY 


. . 

j?r 


i 


f(*) 

B(*) 

A(*) 

P 

W 

Ml 

Y 

PI 

D(*) 


C(*) 


L  (*) 

M(*) 

17,  19,  18 
S 


T 


J,  K,  Kl 


B1 

A1 


Interceptor  array 

Type  B  bomber  array 

Type  A  bomber  array 

Number  of  passes 

Scale  factor  for  type  A  bombers 

Size  of  the  iteration  loop 

Flag  used  to  check  for  new  probabilities  array 
Pass  index 

Probability  of  Detection  and  Conversio’.  for 
type  B  by  pass 

Probability  of  Detection  and  Conversion  for 
type  A  by  pass 

Probability  of  kill  for  type  B  by  pass 

Probability  of  kill  for  type  A  by  pass 

Low,  high,  and  step  size  of  I/B  ratio 

Total  surviving  bombers  of  type  B  (summed  after 
each  iteration) 

Total  surviving  bombers  of  type  A  (summed  after 
each  iteration) 

Loop  indices 

Total  type  B  bombers  in  each  cell 
Total  type  A  bombers  in  each  cell 
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VARIABLE  LIST  FOR  DLMNTY — Continued 

K2,  K3,  X  Dummy  indices  used  in  interceotor  assignment 
Ql  Total  bombers  per  cell 

Q2  Total  bombers  used  in  all  iterations 

all  iterations 
all  iterations 


Q3 

Q4 


Total  type  B  bombers  used  in 
Total  type  A  bombers  used  in 
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Sample  Runs  of  DLMNTY 


3. 


\ 


WHAT  ARE  PASSES#  BOMBER  SCALE  FACTOR#  AND  L00P  SIZE?! 1 .2# 100 
NEW  PR08ABILITIES  ARRAY? CYES-! #N0-O> !  1 

PASS  1  P(DC)#P(K)  (BOMBER  A)#P<DC)#P<K>  (B0M8ER  8> ! .5* .5# *7# .7 

I/B  RANGE  <L0W#HIGH#AND  STEP)?!  .2*2*  *2 

I/B  PROBABILITY  0F  SURVIVAL 


ALL  80MBERS 

TYPE  A 

TYPE  B 

•  20 

.940 

.951 

.918 

•  40 

.878 

.906 

.823 

.60 

.814 

.859 

.724 

•  80 

.764 

.809 

.674 

1  .00 

.728 

.783 

.619 

1  .20 

.667 

.732 

.538 

1  .40 

.625 

.684 

.508 

1  .60 

.61  1 

•  691 

.44° 

1  .80 

.557 

.633 

.405 

2.00 

.525 

.60S 

.366 

WHAT  ARE  PASSES#  BOMBER  SCALE  FACTOR#  AND  L00P  SIZE? ! I #2# 1  DO 
NEW  PROBABILITIES  ARRAY?<YES-1#N0-O)!i 

PASS  1  P<DC)#P<K)  CB0MBER  A)  #P(DC > #P(K)  (30MBER  B ) ! .7 , . 7# . 5# . 5 
I/B  RANGE  <L0W#HIGN#ANO  STEP)?! .2, 2#. 2 

I/B  PROBABILITY  0F  SURVIVAL 


ALL  90MBERS 

TYPE  A 

TYPE  B 

.20 

.923 

.9)0 

.950 

.40 

.848 

.817 

.910 

.60 

.773 

.739 

.840 

.80 

.720 

.673 

.814 

1  .00 

.660 

.598 

.783 

1  .20 

.622 

.562 

.742 

1  .40 

.576 

.512 

.705 

1.60 

.528 

.449 

.686 

1  .80 

.496 

.421 

.645 

2.00 

.455 

.373 

.617 
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WHAT  ARE  PASSES#  B0MBER  SCALE  FACT0R#  AND  L00P  SIZE?! 2#2# 100 
NEW  PR0BABILITIES  ARRAY? C YES -1 #N0-O) ! 1 

PASS  1  P (DC ) #P  CK)  CB0MBER  A)#P <DC >#P CK >  CB0MBER  B > ! . 5# . 5# . 7# . 7 

PASS  2  PCDC)#PCK)  CB0MBER  A)#PCDC >#PCK>  (B0MBER  B > ! .3# *3# . 5# . 5 

I/B  RANGE  <L0W# HIGHLAND  STEF)?! *2#2#.2 

1/8  PROBABILITY  0F  SURVIVAL 


.  :■ 

ALL  BOMBERS 

TYPE  A 

TYPE  B 

>  ■- 

•  20 

•  909 

.937 

•  854 

•  40 

•  814 

.867 

.707 

?  ; 

•  60 

.737 

•  806 

.600 

•  80 

•  66! 

.729 

.  526 

1  .00 

.579 

.663 

.409 

1  .20 

.539 

.651 

.315 

l 

1  .40 

.474 

.574 

.275 

1 

1  *60 

.433 

.549 

.201 

1  .80 

.375 

.494 

.138 

;  | 

2.00 

•  343 

•  459 

.110 

:  5 

WHAT  ARE  PASSES# 

80M8ER  SCALE  FACTOR#  ANO 

L00P  SI2E7IO 

#0#0 

*  The  0,0,  0  input  designates  exit  from  DLMNTY. 
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SAMPLE  RUNS  USXUG  DLMNTP — Continued 
(Merge  DLMNTZ,  DLMNTY,  DLMNTP) 
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IV.  TWO  WAVE  BOMBER-INTERCEPTOR  MODEL— TOOTH 


A.  Purpose 

The  purpose  of  the  TOOTH  series  of  computer  models  is 

to  fill  the  gap  between  simple  Monte  Carlo  models  which 

relate  probabilities  of  l/B  ratios  and  detailed  simulation 

3 

models  involving  basing,  C  ,  geography,  and  operational 
specifics.  The  TOOTH  models  can  be  used  with  basic  inputs 
for  rapid  excursions  or  for  use  with  a  manual  war  game  to 
exercise  the  random  parts. 

B .  Basic  Concept  of  TOOTH 

The  basic  concept  involves  the  penetration  by  four  types 
of  bombers — two  carrying  gravity  weapons  and  two  carrying 
standoff  weapons  (ASMs) .  Defense  is  provided  in  two  waves 
by  two  types  of  interceptors.  For  the  first  attack  wave, 
it  is  assumed  that  there  is  a  degrading  ECM  environment  and 
the  C3  only  provides  a  raid  count  by  cell  number.  The  inter¬ 
ceptors  then  attack  by  committing  a  certain  number  per  cell. 

In  the  second  wave,  it  is  assumed  standoff  weapons  have  been 
launched  and  that  a  bomber  count  has  been  made.  The  weapons 
of  interceptor  type  2  attack  the  ASMs  and  bombers  as  resources 
allow.  A  tally  is  made  of  surviving  bombers  at  the  end  of 
each  wave  and  of  the  number  of  surviving  ASMs  of  each  type 
after  the  second  wave. 

The  program  is  a  stylized  war  game  because  no  geographi¬ 
cal  aspects  are  modeled.  However,  the  operational  aspects  of 
the  game  are  accounted  for  in  the  inputs  to  the  program.  This 
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manual  "homework"  is  used  to  determine  the  number  of  intercep¬ 
tors  that  can  be  brought  to  bear  in  the  first  wave,  the  kind 
of  weapons  they  carry,  and  such  details  as  turn-around.  For 
example,  type  1  interceptors  can  be  added  to  the  second  wave 
"pool"  of  interceptors  to  represent  that  fraction  of  inter¬ 
ceptors  that  could  be  readied  for  the  second  wave  attack. 

Note  that  the  bombers  do  not  release  ASMs  until  after  the 
first  wave  of  interceptors  has  attacked.  This  is  a  rough 
approximation  of  the  geometric/geographic  aspects  of  air 
defense. 

In  this  model,  we  consider  ASMs  to  be  attackable  targets 
and  allocate  interceptors  to  the  bombers  and  ASMs  on  a  random 
basis . 

C .  Program  Description 
1.  Input 

The  following  numbers  are  requested  from  the  user  to 
run  TOOTH: 

•  Number  of  bomber  cells  (Type  1-4),  C(l),  C(2),  C(3), 
and  C (4) . 

•  Number  of  interceptors/cell  and  passes  for  the  1st 
and  2nd  wave  (II,  PI,  12,  P2) ,  where  II  is  the 
number  of  interceptors  of  type  1  committed  to  each 
cell  on  the  1st  wave,  PI  is  the  number  of  passes 
of  type  1  interceptors  against  the  penetrators, 

12  is  the  total  number  of  interceptors  of  type  2 
assigned  to  the  bombers  and  ASM's  for  the  2nd  wave, 
and  P2  is  the  number  of  passes  the  type  2  interceptors 
make  against  the  bombers  and  ASM's. 
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The  following  probabilities  are  requested  in  this  order: 

•  K(  1, 1) ,  K(2, 1)  ,  K( 3, 1)  ,  K(4, 1)  ,  where  I  is  the  index 
for  the  number  of  passes  by  type  1  interceptors  on 
the  1st  wave.  Therefore,  K(3,2)  would  be  the  proba¬ 
bility  of  kill  of  the  type  1  interceptors  against 
type  3  bomber  on  the  2nd  pass. 

•  L  (1, 1)  ,  L(2,I),  L(3,I),  L(4,I),  where  I  is  the  index 
for  the  number  of  passes  by  the  type  2  interceptors 
on  the  2nd  wave. 

•  K3,  K4,  where  K3  is  the  P^  of  the  type  2  interceptor 
against  the  type  1  ASM  and  K4  is  the  of  the  type  2 
interceptor  against  the  type  2  ASM. 

•  R3,  R4,  where  R3  is  the  probability  of  a  successful 
type  1  ASM  launch  and  R4  is  the  probability  of  a 
successful  type  2  ASM  launch. 

The  loop  size  Mi  is  to  be  specified  and  indicates  the  number 
of  times  the  game  is  to  be  played. 

2 .  Simulation 

The  game  is  played  Ml  times  to  improve  statistics;  the 
user  specifies  Ml  each  time.  The  results  are  averages  of 
these  Ml  games. 

In  each  game,  the  bomber  cells  are  chosen  for  attack, 
one  at  a  time.  The  number  of  bombers  per  cell  is  randomly 
chosen  between  1  and  10.  The  number  of  bombers  surviving 
the  first  wave  of  attack  is  accumulated.  A  bomber  is  con¬ 
sidered  killed  or  not  killed  by  a  random  choice  based  on  the 
probability  of  kill  involved.  The  number  of  interceptors 
committed  is  already  known,  so  the  terminal  engagement  i.s 
repeated  until  all  passes  have  been  made  or  all  bombers  are 
killed.  Surviving  bombers  of  each  type  are  tallied  for  the 
second  wave. 
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For  the  second  wave,  it  is  assumed  that  a  count  of  the 
surviving  bombers  of  each  type  is  known.  Also,  the  number 
of  ASMs  launched  is  known.  The  program  allows  those  penetra- 
tors  which  carried  ASMs  to  continue  with  the  raid  for  detec¬ 
tion  and  interception.  The  program  assumes  ASMs  could  be 
attacked  from  the  total  pool  of  type  2  interceptors.  The 
interceptors  are  uniformly  committed  to  bombers  and  ASMs 
for  the  specified  number  of  passes  till  the  supply  is  exhausted 
or  all  bombers  and/or  ASMs  are  killed. 

After  the  specif iod  number  of  passes  on  the  second  wave, 
the  surviving  bombers  and  ASMs  are  accumulated  and  the  entire 
game  is  reinitialized  and  is  repeated  until  Ml  runs  have  been 
made . 

3.  Output 

The  output  totals  are  averages  from  the  Ml  games.  Since 
only  cells  were  specified  on  input,  the  program  ouputs  the 
average  size  of  the  raid  in  terms  of  total  penetrators.  It 
is  assumed  that  the  averages  of  each  bomber  type  are  pro¬ 
portional  to  the  number  of  cells  of  each  type. 

The  summary  of  results  is: 

-1  number  of  bombers  in  attack  and  of  each  type  (1-4) 

-2  number  of  bombers  surviving  1st  wave  of  attack  by 
type  (1-4) 

-3  number  of  bombers  surviving  2nd  wave  of  attack  by 
type  (1-4) 

-4  ASMs  launched  by  bomber  types  1  and  2 

-5  ASMs  surviving  by  types  1  and  2. 
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Coding 


•  Program  List  for  TOOTH 

•  Instructions  for  TOOTH 

•  Variable  List  for  TOOTH 


PROGRAM  LIST  FOR  TOOTH 


100  DIM  ZUO>10>*Y<10*iQ>*AaO#10>*BUQ#10> 

1000  PRINT  "NUMBER  0F  B0M3ER  CELLStTYPES  1  T0  4)?"! 

1010  INPUT  C<1 >*CC2)*C(3>*C<4> 

1020  PRINT  "NUMBER  INTERC/CELL  AND  PASSES*  1ST  AND  2ND  WAVE?"! 
1030  INPUT  1 1 *P1 * I2*P2 

1040  PRINT  "NEW  PROBABILITIES  ARRAY  CYES- 1 *N0-2 > ?"! 

1050  INPUT  Y 

1060  IF  Y*2  THEN  1300 

1070  FOR  1=1  TO  PI 

1080  PRINT  "PASS"I J"1ST  WAVE  P(K)  VS  EACH  BOMBER  TYPE?"! 

1090  INPUT  K(1*I>*K(2*I)*K(3«I)*K£4*I) 

1100  NEXT  I 
1110  FOR  I«1  TO  P2 

1120  PRINT  "PASS" I !"2ND  WAVE  P<K>  VS  EACH  BOMBER  TYPE?"! 

1130  INPUT  LU*I)*L<2*I>*LC3*I>*L<4*I) 

1140  NEXT  I 

1150  PRINT  "P (K )  VS  ASM-S* 1ST  TYPE*2ND  TYPE?”! 

1160  INPUT  K3*K4 

1190  PRINT  "ASM  LAUNCH  RELIBILITY  -1ST  TYPE*2N0  TYPE?”! 

12U0  INPUT  R3*R4 

1300  PRINT  "LOOP  SIZE?"! 

1310  INPUT  Ml 

1500  MAT  B=ZER 

1600  MAT  A=ZER 

1800  FOR  M* 1  TO  Ml 

2060  NAT  Z=ZER 

2070  MAT  Y*ZER 

2100  FOR  T=1  TO  4 

2105  IF  C(T)=0  THEN  2390 

2110  FOR  1=1  TO  C<T> 

2120  B9=INTOO*RND(Om> 

2130  Z<T*3)=Z<T#3)*B9 

2135  IF  11=0  THEN  2370 

2140  15=11 

2200  FOR  P=1  TO  PI 

2210  IF  B9>»I 1  THEN  2300 

2220  I 5=B9 

23Q0  FOR  L= 1  TO  15 

2310  IF  RND(0)>KCT*P)  THEN  2350 

2320  B9=B9-1 

2330  IF  B9=0  THEN  2370 

2350  NEXT  L 

2360  NEXT  P 

2370  ZCT»1 )=Z<T*1 )+B9 

2380  NEXT  I 

2390  NEXT  T 
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PROGRAM  LIST  FOR  TOOTH — Continued 


2500  Y<t*3)=INTCR3*Z(l*l)> 

2510  YC2*3>*INT<R4#Z(2*l >> 

2520  YU*n»Y<l#3> 

2530  Y<2*1 >*Y<2*3) 

2600  F0R  J*1  T0  4 
2610  Z<J*25*Z<J*I> 

2700  NEXT  J 

2900  A9«Y<1 *1)+Y<2*1)+Z(1*2)'*'Z(2*2)*Z<3*2>  +  Z<4*2) 


3000  FOR  1*1  T0  12 
3010  9*INTC6*RND<0)+1 ) 

3020  G0T0  3090*3190*3290*3390*3490*3590  0N  Q 

3090  IF  Y< 1 « 1 >»0  THEN  3010 

3100  IF  RN0(Q)>K3  THEN  3700 

3110  YC1,1 >*Y<1*1>-1 

3120  A9*A9-1 

3130  GOTO  3700 

3190  IF  YC2* i )«0  THEN  3010 

3200  IF  RND <0  )  >K4  THEN  3700 

3210  Y(2*l)aY(2*l)-l 

3220  A9=A9-1 

3230  GOTO  3700 

3290  IF  ZC1 *2)»0  THEN  3010 

3300  IF  RN0(0)>L(1*S)  THEN  3700 

3310  Z  <  1  *2  )»Z< 1 *2 ) - 1 

3320  A9*A9-1 

3330  GOTO  3700 

3390  IF  Z  <2*2  )»0  THEN  3010 

3400  IF  RN0<0)»L<2*1 >  THEN  3700 

3410  Z<2»2)bZ(2*2>-1 

3420  A9-A9-1 

3430  GOTO  3700 

3490  3F  Z<3*2)«0  THEN  3010 

3500  IF  RNQ<0 >>L  <3* 1 )  THEN  3700 

3510  ZS3»2>s»Z<3*2)-l 


SS20  A9*A9-S 

3530  G0TS  3700 

3590  IF  Zf4*2)*0  THEN  3010 

3600  IF  AND<0>»lf4#l.)  THEN  3700 

3610  Z<4*2)*Z<4»2>-1 

3620  A9«A9^1 

3630  GOT®  3700 

3700  IF  A9®0  THEN  3800 

3710  NEXT  I 

3800  FOR  T«1  T®  4 

3810  atTs2>*B<T*2>*Z<T*2> 

3820  B<T*3)*B<T#3!'+2<T*3l 
3830  9<T*t  )»8<T*n^Z<t*n 
3840  NEXT  T 
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PROGRAM  LIST  FOR  TOOTH— Ccntinued 


3850  A(1,3)=A(|,3)+Y(!,3) 

3860  A<2,3>=A<2*3)+Y<2#3> 

3870  ACl*l>sA(l,l )+Y(l>l ) 

3880  A(2«1)=A<2«1 )+Y<2*l ) 

3900  NEXT  M 

4000  8*CB(l#3)+8<2jT)+B<3i>3>+B(4*3>>/MI 
4010  A*(A<1»3)+A<2>3))/Mi 

SOOO  PRINT  “TOTAL  ATTACK  SIZE  BOMBERS'*! 8l"ASM-S“l A 
5100  PRINT  “ATTACK  SIZE-BY  BOMBER  TYPE'* 

5110  PRINT  BU*3>/MlJBC2«3)/Mne<3*3>/MUB«4*3>/Ml 
5200  PRINT  “BOMBERS  SURVIVING  FIRST  WAVE1* 

5210  PRINT  3(!*l)/MlJB(2»l)/MlJBC3«l)/M!}B(4ifl)/Mt 
5300  PRINT  "BOMBERS  SURVIVING  SECOND  WAVE" 

5310  PRINT  B<1,2)/M11B(2,2>/M1JB<3#2)/MHB<4*2>/M1 
5400  PRINT  "ASM-S  LAUNCHED-BY  BOMBER  TYPE" 

5410  PRINT  A<t  #3)/MllA(2i*3>/Mi 
5500  PRINT  "ASM-S  SURVIVING" 

5510  PRINT  AU  »\ )/MllA(2#J  )/Mt 

7000  PRINT  “END  OF  GAME?  <YES~l)“l 

7010  INPUT  G 

7020  IF  G*1  THEN  9999 

7100  GOTO  1000 

9999  END 


INSTRUCTIONS  FOR  TOOTH 


1000  -  1310  Input  Section. 

1000  -  1010  The  number  of  bomber  cells  for  each  type 
of  interceptor.  NOTE:  Types  1  and  2  are 
ASM  carriers. 

1020  -  1030  II  is  the  number  of  type  1  interceptors  per 
cell  committed  in  the  final  wave  and  PI  is 
the  number  of  passes  made.  12  is  the  total 
number  of  type  2  interceptors  in  the  second 
wave,  since  there  are  no  "cells"  in  the 
second  wave.  P2  are  the  number  of  passes 
made  by  the  type  2  interceptor. 

1040  -  1200  If  a  new  array  is  to  be  set  up,  then  the 

request  is  made  for  the  probability  of  kill 
of  each  interceptor  against  each  bomber  type 
for  each  pass.  Also,  the  program  requests 
the  P  and  reliability  of  launch  of  each 
ASM  tfpe. 

1300  -  1310  Number  of  iterations. 

1500  -  1600  Bomber  arrays  are  initialized  to  zero. 

1800  -  3900  Outer  iteration  loop. 

2060  -  2070  Bomber  arrays  2eroed  for  each  iteration . 

2100  -  2390  Raid  Attack  Loop  (by  type). 

2105  If  no  cells  of  type  T,  next  cell. 

2110  -  2380  Cell  Attack  Loop  (by  #  colls) , 

2120  Pick  number  1  through  10  randomly  for 
numbers  of  bombers  in  the  cell. 

2130  Sum  Bombers  by  type  within  iteration. 

2135  If  no  interceptors/cell,  skip  attack. 

2200  Pass  Loop. 
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INSTRUCTIONS  FOR  TOOTH — Continued 


2210  If  more  interceptors  than  bombers 
set  number  of  interceptors  to  num 
ber  of  bombers  left  (2220) . 

2300  Interceptor  Attack  Loop. 

2  310  Check  for  Kill. 

2320  If  kill,  decrease. 

2330  Check  for  all  bombers  killed  in 
current  cell. 

2350  Next  Interceptor. 

2360  Next  Pass. 

2370  Sum  bombers  surviving  first  wave. 

2380  Next  Cell. 

2390  Next  Bomber  Type. 

2500  -  2510  Total  ASM's  launched  (based  on  surviving 
bombers  and  launch  reliability) . 

2520  -  2530  Initialize  ASM's  for  second  wave. 

2600  -  2700  Initialize  bombers  for  second  wave. 

2900  Total  targets  for  second  wave. 

3000  Interceptor  attack  loop. 

3010  Pick  random  number  1  to  6 . 


3020  Branch  to  subset  of  instructions  depending 
on  targets. 


3090  -  3130 
3190  -  3230 


Attacks  on  the  two  types  of 
ASMS. 


3290  -  3330 
3390  -  3430 
3490  -  3530 


Attacks  on  the  four  types  of 
bombers.  In  each  case,  kill 
check  is  made,  sum  is  decreased, 
and  exit  to  3700. 


3590  -  3630 


3700  If  all  killed,  go  to  summarizing 
section  (3800)  . 
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INSTRUCTIONS  FOR  TOOTH— Continued 


3710 

Next  interceptor. 

3800 

-  3840  Total  Bombers  and  those  surviving  each 
wave  are  added  to  the  total  for  all 
iterations . 

3850 

-  3880  Total  ASM's  launched  and  those  surviving 

the  2nd  wave  are  totaled  for  all  iterations 

3900 

Next 

Iteration. 

4000 

-  4010 

Bombers  and  ASM's  averaged. 

5000 

-  5510 

Print  out  section. 

7000 

-  7100 

Continuation  request. 

9999 

END 

h 

’  i; 
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VARIABLE  LIST  FOR  TOOTH 


I  * 


:  1  .-: 
y  i 


Z  (**)  Number  of  bombers  for  each  iteration 
B  (**)  Running  total  of  bombers 
Y  (**)  Number  of  ASM's  for.  each  iteration 
A  (**)  Running  total  of  ASM's 


NOTE:  In  the  above  arrays,  the  first  position  represents 
the  type  of  bomber  or  ASM.  The  second  position  designates 
the  situation  as  follows: 

-1  those  which  survive  the  first  pass 
-2  those  which  survive  the  second  pass 
-3  total  at  the  start  of  the  game  . 

Thus: 

Z  (3,2)  is  the  number  of  bombers  of  type  3  which 

survive  the  second  pas.5?  for  that  iteration. 

A  (1,3)  is  the  total  number  of  ASM's  of  type  1 

which  have  been  launched  in  all  iterations 
up  to  that  point. 

The  A' s  and  B'a  are  the  running  totals  of  the  individual  Z's 
and  Y'a.  At  the  end  of  all  iterations,  they  are  averaged  to 
find  the  results  of  the  "average"  game. 

C  (*)  Number  of  cells  of  each  type  (1-4) 

11  Number  of  irterceptorn/cell  committed  first  wave 

12  Total  number  of  interceptors  second  wave 
PI  Number  parses  first  wave  interceptors 

P2  Number  passes  second  wave  interceptors 

K(**)  Kill  probability  array  for  first  wave?  1st  position  is 
bomber  type,  2nd  position  is  pass. 
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VARIABLE  LIST  FOR  TOOTH— Continued 


mhwww 


I 

i 


L  (**)  Same  as  K  (**)  for  second  wave 
K3  P^  against  1st  ASM  type 

K4  P  against  2nd  ASM  type 

K 

R3  Reliability  of  launch  of  1st  ASM  type 

R4  Reliability  of  launch  of  2nd  ASM  type 

Ml  LOOP  size  for  iteration 

T  Bomber  type  ■y 

B9  Number  bombers  within  the  current  cell 

A9  Total  targets  for  second  wave 
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5.  Sample  Runs  of  TOOTH 


NUMBER  0F  B0MBER  CELLSCTYPES  1  T0  4)?? 
NUMBER  INTERC /CELL  AND  PASSES*  1ST  AND 
NEW  PROBABILITIES  ARRAY  < YES- 1 *NQ-2 > ? ? 


PASS  1  1ST  WAVE  P (K )  VS 
PASS  2  1ST  WAVE  PCK)  VS 
PASS  3  1ST  WAVE  PCK)  VS 
PASS  1  2ND  WAVE  PCK)  VS 
PASS  2  2ND  WAVE  PCK)  VS 
PASS  3  2ND  WAVE  PCK)  VS 
PCK)  VS  4SM»S*1ST  TYPE*2ND 
ASM  LAUNCH  RELIBILITY  -1ST 


EACH  BOMBER 
EACH  BOMBER 
EACH  BOMBER 
EACH  BOMBER 
EACH  BOMBER 
EACH  BOMBER 
TYPE??  .3**1 
TYPE.2ND  TYPE?? 


15*10* I  5* ID 

2ND  WAVE??  3*3*300*3 

I 

TYPE??  .7*.S*.7*.5 
TYPE??  *6»*4**6**4 
TYPE??  .6*. 4*. 6*0 
TYPE??  .6*. 4*. 6*. 4 
TYPE??  .5*. 3*. 5*. 3 
TYPE??  .5*. 3*. 5*0 


•  8*  •  8 


LOOP  SIZE??  100 

TOTAL  ATTACK  SIZE  BOMBERS  268.14  ASM-S  33.61 
ATTACK  SIZE-BY  BOMBER  TYPE 
81.29  53.79  78.91  54.15 

BOMBERS  SURVIVING  FIRST  WAVE 
19.6  23,47  18.08  31.31 

BOMBERS  SURVIVING  SECOND  WAVE 
8.00000E-0?  2.88  .18  6.22 

ASM-S  LAUNCHED-BY  BOMBER  TYPE 

15. 27  18.34 
ASM-S  SURVIVING 

1.27  11.37 

END  OF  GAME?  CYES-D?  0 

NUMBER  ”'0?  ^OmAIPhT’C  eiLS  <  TYPES  1  TO  4  >  ? '?  I  5  #  10*15*10  ~ 

NUMBER  1NTERC/CELL  AND  PASSES*  1ST  AND  2ND  WAVE??  3*1*100*3 
NEW  PROBABILITIES  ARRAY  CYES- 1 *N0-2 ) ??  2 


LOOP  SIZE??  100 

TOTAL  ATTACK  SIZE  BOMBERS  274,9?  ASM-S  76.15 
ATTACK  SIZE-BY  BOMBER  TYPE 
82.82  55.17  83.44  53.49 

BOMBERS  SURVIVING  FIRST  WAVE 


54.37  41.74  54.88  39.9 

BOMBERS  SURVIVING  SECOND  WAVE 
44. n  35.07  44.81  33.22 

ASM-S  LAUNCHED -BY  BOMdER  TYPE 
43.14  33.01 

ASM-S  SURVIVING 


37.86  31.11 

END  OF  SAME?  CYES-l)?  0 
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NUMBER  0E  BUNDER  CEU-S6TYRES  1  TO  4)??  1  5.  1 0*  1  S,  1 0 
NUBBER  INTERC/CFLL  AND  PASSES*  1SI  AND  2ND  WAVE??  3*1*100*1 
NEE  RRUfcAbILl  TIES  AhhAY  <  Y  b  S-  1*  Nk»-2  >  (  ?  2 
LUOR  SIZE??  100 

TOTAL  ATTACK  SIZE  hUMBb.KS  275.36  AM*.-?  75.31 

ATTACK  SIZE- BY  DUMBER  TYRE 
82.61  54*67  83.56  54.52 

BOBBERS  SURVIVING  HIRST  WAVE 
54*19  40.94  54.97  41.05 

BUNDERS  surviving  second  wave 

44.14  34*27  45*31  34 

ASM-S  l.AUNCRED-bY  BOBBER  TYRE 
42.95  32.36 

ASB-S  SURVIVING 
37.82  30.64 

EM)  UK  DANE?  (YES-1)?  0 _ _ _ _ 

NUMBER  OH  bUMbbR  CELLSt TYRES  1  TJ  4)??  15*10.1^*10 

NUMBER  1  iv'i  ERG/C  F.Ll  AND  RaSSE  S*  1ST  AND  Pi\IV  WAVF??  1*3*  |i "0*3 

NEW  RRUDaBILITI  ES  ARRAY  (YES-  1  *  NO- 2) ?  ?  P 

LUuR  SIZE??  100 

TUT/4.  ATT  Aon  MZb  BunDERS  279.06  AS.V-S  7R.7 
ATTACK  SIZE'-tY  DUMBER  TYRE 
b>2  •  6 1  55.91  84.57  55.97 

bUt*. E'ER?  SURVIVING  HIRST  WAVE 
55.42  43.8  4  57.48  47.26 

DOMbbRS  SURVIVING  SECOND  WAVE 
45*66  3  i  .  1 9  47.14  4 1  .  1 3 

aSm-R  LAUNCHbU-BY  f.AvhER  TYRE 
44. 01  34.6V 

ASM-S  SURVIVING 
39.04  32.99 

END  OH  GAMH?  (YtS-P?  C 

nUM4R'  Ub  UoBUfch  CELUSn VTTs  ’ri j  '47  >?  IS.ic.  l  s'io 
NUMtitR  InTERC/GELL  AM:  PASSES*  1ST  AND  PM-  WAVE??  1*1*100,1 
NEW  RRUbAbiLi  1 1 ES  ARRAY  < Y b S- 1  * NJ-p 7 ? t  ? 
lUOR  SIZE??  1UC 

TOTAL  All  AON  SIZE  BONDERS  274.54  A  SB  -  ?  96.0* 

ATTACK  SIZE -bY  DUGGER  TYRh 
8  2*2.2  34«lt>  h3»13  55»G1 

RUBBERS  SURVIVING  HIRST  WAVE 
7  1  .  oV  49*2  7  7  2.67  4‘> .  R 

BONDERS  SURV 1  U  1  nG  Sfc  GUN U  WaVH 
61. £6  42.55  6P.66  42*84 

ASB-S  LMUNCmHj-RY  bob  i'ER  mb 
56*99  39.0V 

aSk  —  S  SURVIVING 
52* t 5  37.4 

END  UH  U AH. t i  (YES-1)?  1 
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V.  EPILOGUE 

The  three  models  described  in  this  Note  have  been  in 
use  at  ANSER  for  a  number  of  years  in  many  projects.  Their 
results  have  proved  to  be  comparable  with  those  of  large- 
scale  simulations  and  expected -value  models,  even  though 
the  large-scale  simulations  show  much  more  detail  and  the 
expected-value  models  less  detail. 

The  three  models  are  economical.  Typical  runs  cost 
only  a  few  dollars.  Many  excursions  have  been  run  for 
comparisons  of  interceptor  missiles  and  fire  control  systems. 

The  keys  to  the  results  lie  in  the  homework.  Essentially, 
these  programs  relieve  the  analyst  of  the  bookkeeping  associated 
with  a  small  war  game  and  give  him  quick  feedback  on  varia¬ 
tions  of  the  input  parameters. 

There  have  been  many  modifications.  The  models,  being 
simple  and  short,  allow  room  for  imagination  and  variation. 
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A  A  C  T 


This  Note  discusses  three  time-sharing  computer  programs  which  are 
quick  tools  for  air  defense  engagement  analysis.  The  programs  are 
based  on  a  Monte  Carlo  method.  The  model  MONTYX  evaluates  air  defense 
interceptor  configurations  by  determining  their  effectiveness  against 
a  bomber  cell  in  the  terminal  phase  of  an  intercept.  The  model  DLMNTY 
expands  the  capability  of  MONTYX  to  a  situation  where  two  bomber  types 
are  contained  within  the  bomber  cell.  The  results  are  bomber  probabil 
ities  of  survival  for  various  interceptor  to  bomber  ratios.  The  model 
TOOTH  involves  the  penetration  by  four  types  of  bombers--two  carrying 
gravity  weapons  and  two  carrying  standoff  weapons.  The  defense  is 
provided  by  two  types  of  interceptors  in  two  waves.  In  the  first  wave 
only  raid  count  by  cell  is  provided.  The  second  wave  assumes  that  a 
bomber  count  has  been  made  and  standoff  weapons  have  been  launched. 

The  results  indicate  the  numbers  of  bombers  of  each  type  surviving 
each  wave  of  defense  as  well  as  the  number  of  standoff  weapons  launched 
and  those  which  survive  the  second  wave  of  defense. 
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